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The crucial formulas to determine the non-chaotic states in the rf-biased non- 
linear oscillators are derived from the numerical experiments. The nature of these 
formulas, which depends on symmetrical properties of the potential well, in terms 
of the driven-frequency as a function of the damping constant k is investigated. All 
these ones provide crucial guide posts to check which kinds of solutions (simple or 
complicated) can be tailored in the dissipative rf-biased nonlinear oscillators, respec- 
tively. 
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1. INTRODUCTION 



In the nonlinear oscillating systems with an external driven force, the equation of motion 
can be expressed as a general form: 



where k denotes the damping constant, F and uj represent the amplitude and the fre- 
quency of the external periodic driving force respectively, and f(x) is the nonlinear term. If 
f(x) can be expressed as x + + ■ ■ the Eq.(l) becomes as the dynamic motion of Duffing 
oscillators. Whereas if f{x) — sinx, the equation is called the resistant-shunted Josephson 
junction model. Following previous work^~^, the complicated solutions to the Eq.(l), includ- 
ing period-2^, chaotic, period-3 x 2^ and period-m x 2^ ones {k is a positive integer and m a 
prime), have been formed in the parameter region where cu and k are rather small. Moreover, 
the non-chaotic solutions to the Eq.(l), including symmetric and asymmetric period- 1 ones, 
dominate in a manner of large k and large lu. In this case, no matter what F is, there is only 
symmetry-breaking bifurcation observed. In fact, an occurrence of non-chaotic states was 
first indicated by McDonald and Plischke^ in the study of the resistant-shunted Josephson 
junction (RSJ model^°'^^) with the following criterion: 



i.e., the damping constant is much larger than a unit. And then subsequently Kautz 
pointed out that chaos does not occur for 



According to their suggestions, there are only simple or stable states survivable in some 
confined region of the parameter space of the higher damping and/or the faster driving 



X + kx + f{x) — F sin cut 



(1) 



1 « e 



(2) 



u » k 



(3) 
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frequency cases. Although these criteria can enhghten the parameter regions of the comph- 
cated motion quahtatively, a quantitative resolution of this problem need explore yet. In 
other words, the results described above immediately raises an intriguing problem : can one 
exactly predict where or how a nonlinear oscillating system will display non-chaotic states 
in terms of chosen parameters cu and k quantitatively. 

To this end, herein, we will attempt to deduce certain crucial formulas in terms of the 
parameters, cu and k, to quantitatively distinguish the existence of chaos and the occurrence 
of non-chaotic states in Eq.(l) not only in the RSJ model but also in Duffing oscillators. In 
what follows the method of derivation of the formulas with the aid of the phase diagrams, 
bifurcation diagrams, state diagram and so on will be described. As compared with the RSJ 
model and the generalized Duffing oscillators, the physical implications of the symmetrical 
properties of potential wells related to the character of the formulas will be examined. 
Concluding remark will be reported as well. 

2. The Numerical Simulation 

As to derive the crucial formulas to distinguish the existence of chaos from the occurrence 
of non-chaotic states in the more general cases of the dissipative nonlinear oscillator, two 
typical nonlinear oscillating systems, RSJ model and Duffing oscillators, are considered at 
present. In the numerical simulation, the fourth order Runge-Kutta algorithm^^ is employed 
to do integration. Based on the integration, the Newton- Raphson method^^ is utilized to find 
the stable and unstable fixed points on the Poincare sections. Two scanning procedures, 
by varying driving frequency at a fixed driving amplitude (a; scanning) and by varying 
driving amplitude at a fixed driving frequency {F scanning), are employed with damping 
constant A; as a fixed parameter to get as much information as possible. Consequently, 
various characteristics of all the thresholds of possible transitions in terms of the controlled 
parameters can be obtained and depicted in a state diagram. 

In each state diagram of the respective systems, for example as shown in Fig.l, there may 
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exist a threshold frequency, a;T, at a fixed damping constant, k. We found that if the driving 
frequency is larger than the frequency, a;^, in this damping case, no matter how large the 
driving force amplitude is the chaotic solutions cannot exist. Then, as varying the damping 
constant k to the other fixed value we can obtain another ujt with the same process. Finally, 
to repeat the process again and again, the relations between uit and k can be established 
in this way after all. These relations allow us to derive a border line and then to formulate 
a closed mathematical form in the uj — k space, by which the complicated solution region 
from the simple stable period- 1 solution one can be differentiated. The detail contents of 
such UJ — k relations of the formulas for the RSJ model and the Duffing oscillators will be 
discussed below. 

A. RSJ Model 

The rf-biased RSJ model is given as: 

X -\- kx -\- sin a; = Fsina;^ (4) 

where k is the damping constant, F is the amplitude of rf driving force and uj is the 
frequency of rf driving force. Eq.(4) describes the equation of motion of the phase difference 
of the junction. The equation can also be visualized as the dynamic motion of the swept 
angle of pendulum^'^^'^^, charge density wave^, a parametric amplifiers^^'^^, et. al.. By 
means of the numerical method described above, the features of the state diagram related 
to various solutions of the states^^^^ in Eq.(4) will be presented as follows. For the sake of 
illustration, only seven branches (shaded regions) are denoted in the F — uj states diagram as 
A; = 1.0 as shown in Fig.l. There are two kinds of solution regions. One contains the simple 
solution including the period- 1 oscillatory (symmetric and symmetry-broken) and the phase- 
locked traveling ones. And the other contains the complicated solution including period-2*, 
chaotic, period-3 x 2^ , period-m x 2^ (k is a positive integer and m a prime) and the phase- 
unlocked traveling ones. The shaded region of the branches indicate the existence of the 
complicated solutions. Outside these regions, especially in the region where the frequency is 

4 



above ut, there does exist the simple solutions without observing the chaotic ones. In fact, 
the scenario of a sequence of bifurcation in each branch along the F-axis at a fixed frequency 
cu, is mostly the same, as shown in Fig. 2. They mainly contain period- doubling cascades, 
chaos, intermittency, crisis, and reverse period-doubling cascades. The period-3 windows 
and sometimes the period-5 windows very often burst out within the chaotic events. Since 
the similar scenarios of bifurcation appear itself again and again in each branch, it implies 
the self-similarity of the solutions embedded in the RSJ model clearly. The greater details 
of such self-similar properties will be seen subsequently. 

As also shown in Fig.l, the territory of the each branch is bounded by the terrain of 
the symmetry-breaking (SB) solution, respectively. There are the regions of the symmetric 
solution existed between the regions of the SB solutions as well. Moreover, both the fre- 
quency of symmetric and SB solutions are period- 1. The shaded regions of the branches will 
shrink toward smaller and smaller, closer and closer and then finally merge together when 
cu is decreasing and F is increasing further and further. We have find that in the left-upper 
corner of the Fig.l, actually, there are still much more branches existed. 

In addition to the evidences of the above observations, the self-similar feature of the 
bifurcation in each branch in the RSJ model can be further illuminated as follows. Firstly, 
as shown in Fig.l, the terrains of each simple solution separated by the branches of the 
complicated regions are marked on a series of consecutive integers respectively from down 
to up. We note that the integers just reveal the corresponding number of the local wells of 
the infinite sinusoidal potential within which the oscillating motion of the simple solution 
executes in the RSJ model. For example, in the case of a; = 0.313, F — 1.20 and k — 1.0, 
the response is oscillating in the n — 2 regime as shown in Fig.l, and corresponds to the 
motion swinging over two local wells of the sinusoidal potential within the interval of one 
period of the rf-biasing as shown in Fig. 3b. Similarly, the phase portraits, as shown in Fig.3, 
illustrate these behaviors apparently from n = 1 to n = 7, respectively, when F is increasing 
from 1.80 to 7.0 for u = 0.313. The situation is further demonstrated from Fig.3g for n = 51 
as well. In the complicated region, the common effect of the response swinging extensively 
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from "n-th" local well to "n + l-th" one results from two stable solutions merging with an 
unstable one which located on the top of the sinusoidal well. In other words, the respective 
asymmetrical solutions in the adjacent wells are merging together into one solution. It 
leads to the swinging extensively one more well after the collision between the asymmetrical 
solutions and the unstable one. 

It is also worthy noting that all the branches in the state diagram at a fixed damping 
constant k are located within the region where the driving frequency is less than the crit- 
ical threshold frequency cot- If the driving frequency is larger than the critical threshold 
frequency ujt, no matter how large the driving force amplitude is, the complicated solutions 
cannot exist. As mentioned before, for a fixed damping constant k, we obtain a threshold 
frequency cut- Consequently, we find quantitatively that the critical threshold frequency ujt 
at the first branch, n = 1 , is dependent on the damping constant k with a simple relation: 

2" 

(5) 

where ujq = 1.450 and ko = 1.681. In the Fig.4a, it shows the relation between the 
critical threshold frequency, ujt, and damping constant, k (both the numerical data, • , and 

the fitting function, ). It gives the clear-cut borderhne to differentiate the regions of the 

complicated solutions and the simple ones. According to this relation, it easily reveals that 
chaos does not occur not only in the overdamped case, k » 1, but also in the underdamped 
case to » k. Although these results have been ever qualitatively suggested by McDonald 
and Plischke and by Kautz previously, we herein can quantitatively predict it from Eq.(5) 
that if the damping constant, k, excesses a threshold value, /cq, then no chaotic behavior 
can be found. Therefore, we emphasize that the Eq.(5) provides a valuable information of 
determination of the non-chaotic states in the RSJ model quantitatively. 

According to the self-similar feature mentioned before, we further conjecture that each 
maximum driving threshold frequency of the n-th branches, uJTm depends on the damping 
constant k as the critical threshold frequency ujt (for the branch n = 1) does. By Eq.(5), 
we guess that it might be scaled as the following: 
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where 5 is a parameter dependent on the damping constant k. In what follows the 
equation will be verified tentatively in order to make it appear plausible. 

Firstly, the respective threshold frequencies of the n-th branches, cotuS, are traced out 
under one fixed damping constant k by the driving amplitude scanning method {F scan- 
ning). Then, changing the damping constant k to another one, a new set of the threshold 
frequencies, cutuS, therefore, can also be found out in the same fashion. After collecting 
all of the UJTn ~ ^ relations depending on various values of the damping constant k, we 
summarize them in Fig.4b. In the log-log coordinate, we find the linear relations between 
the logarithm of the threshold frequency cutu and of the branch number n accordingly. The 
property of ujtti and n accordingly enables us to determine its functional equation : 



(7) 



where the function Uc{k) is dependent on the damping constant k and the exponent d 
is the slope of the straight lines in Fig.4b for the corresponding damping constant k. Note 
that such simple power law relation apparently reveals the nature of renormalization, i.e. it 
reproduces itself upon rescaling. 

In order to verify the existence of Eq.6, we are still necessary to examine the properties 
of u!c{k) in greater details. It immediately raises a question that whether the function cudk) 
for arbitrary n can be taken as the functional form 1 — (k/ko) for n = 1, in 

Eq.(6). To this end, we depict the Fig.4c. The content of the figure is clearly shown in a 
manner that the property of the function cUc (as expressed in •) is consistent with that of 

the functional form UJ'j' — UJq ' 1 — {k/ko) (as expressed in ) successfully. Note, herein, 

that the exactness of Eq.6 is totally a consequence of the fact concerning the emergence of 
the self-similar feature under the nonlinear dynamics through the RSJ model. 

Following the spirit of the physical implication of the renormalization method, we at- 
tempt to explore the features of the exponent 5 to the extent at present. Since the value of 



the exponent S depends on the damping constant k (as in Fig.4d), we find the three different 
cases in the following. For the case of the small damping constant k (< 0.1), the exponent 
S reaches a constant (~2.589) and is nearly independent of the damping constant k. As the 
case of the damping constant k getting larger and larger (>0.5), the increasing of exponent 
5 behaves as hyper tangent form and then reaches to its maximum. Finally, the case of 
right after the maximum, the exponent 5 drops rapidly to zero as the value of the damping 
constant k approaching to the threshold value, ko. All these results of the exponent 5 can 
be utilized to realize the self-similar behavior quantitatively. 



B. Duffing Oscillators 

The generalized rf-biased Duffing oscillator governed by the equation is given as: 

x + kx-{ ^^Fsmut (8) 

ax 

where the overdot denotes the derivative with respect to time t, k is the damping factor, 
and V{x) is an anharmonic potential function. This equation has been utilized to model 
a wide variety of physical systems such as the optical bistability in the multiple-photon 
absorption process, soft and hard springs, buckled beam, four wave interaction, and plasma 
oscillation^^"^^, et. al.. In general, the potential function V{x) is described by : 

V{x) = ^x' + ^x' + Y (9) 

with a, and 7 coefficients. Actually, by means of some transformations^^, it embraces 
four fundamental types of potential: Vi — \x'^ — V2 — + V2, — \x'^ — ^x^, and 
V4 = ^x'^ + jx'^. The features of the state diagrams of these four types are hsted, respectively, 
as follows; 

• For the case of Vi potential {Vi = \x'^ — i^^)^'''^^: 

The state diagram is shown in Fig. 5a. The transition boundaries include hysteresis, 
period-doubling (PD), crisis, and intermittency. In the state diagram, we also note that 



the curve of PD folds back at a frequency ut- The comphcated solutions exist only with 
the parameters w below the threshold value uit- According to our experiments as same as 
that of the RSJ model, the threshold frequency ujt is found to be a function of the damping 
constant k with a simple form, 



where ujq = 2.226, ko = 2.200, ci = 0.101, C2 = 4.011, cg = 4.495, and a = 1.585 (see 
Fig. 5b). ko is required to satisfy uJT{ko) — 0. Moreover, the non-trivial coefficients of odd 
order of (k/ko), c\ and C3, are induced by the asymmetrical property of potential well with 
respect to the origin, while they are equal to zero in the case of symmetrical ones. 

• For the case of V2 potential (1^2 = —^^^ + |a;^)^^'^°: 

Fig. 6a shows the corresponding state diagram. The potential well V2 has two local wells 
separated by a bump in the middle of V2. As regards allowing the motion in one of the 
local wells, the shape of the transition boundaries looks like a swallow tail. The curves of 
period-doubling arc folded. The threshold frequency lot is found to be dependent on the 
damping constant k. And this relation can be fitted to satisfy the following equation: 



where the constants a;o — 3.218 and k^ = 2.050 (see Fig. 6b). If the driving frequency 
is beyond the threshold frequency ujt and the driving force amplitude does not exceed the 
boundary H„p, the complicated solution cannot occur. Eq.ll provides a crucial condition 
to determine whether the solution is simple or complicated. All these observed features in 
this system are the same as those of the RSJ model. With further increasing the excitation 
amplitude up to curve the solution becomes stable with the swing throughout two 
valleys. In this situation the dynamics of the swing closely resembles the case of the infinitely 
bounded potential V4, while the small effect of the bump in the well is negligible. 




(10) 



oJrik) — ujq 1 



k_ 
ko 



(11) 



• For the case of Vs potential {V3 — — |x^) 
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The state diagram for the transition boundaries with the swallow-tailed form is shown in 
Fig. 7a with damping constant k = 0.1. the transitions include hysteresis, symmetry breaking 
(SB), PD, crisis and intermittency. Due to the fact that the potential barrier is finite, the 
response will escape from the well after a sequence of bifurcation of Feigenbaum period- 
doubling and/or intermittency routes to chaos. The threshold frequency, ujt, at which the 
curve of PD folds back is well fitted in term of the damping constant k as the following 
relation: 





' k' 
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(12) 



where the constants loq — 1.528 and ko — 1.985; (see Fig. 7b). With the driving frequency 
higher than the threshold frequency, ujt, only the simple period-1 solutions are obtainable. 



• For the case of V4 potential (V4 = ^x'^ + jx^)^"^ 



The potential is symmetrical, without an inflection point, and is infinitely bounded as 
\x\ 00. The state diagram is shown in Fig.8. The shapes of the transition boundaries are 
no longer like the swallow tail and can be classified into two groups with characteristic shapes 
associated with odd and even resonances. The chaotic solutions exist in the even resonant 
regions (marked as "C" in Fig.8). There is no such threshold frequency cut observed. Also, 
there is no appearance of the similar cut — k functional relation function in a closed form. 



3. Conclusion 

In this work, we have systematically examined the rf-biased Josephson junction and the 
Duffing oscillators in a wide range of the parameter space. The crucial formulas of the 
functional relation, cut vs. k, have been derived to give the clear-cut criteria for identi- 
fication of the existence of the non-chaotic states. Together with our preceding report^^, 
it is worthwhile mentioning that the functional cut — k relation can be obtained in those 
potentials which possess the inflection points only. Furthermore, for those potential wells 
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possessed symmetrical property with respect to the origin of the potential, the a;^ — A; re- 
lation is in the form of 
the form of ujt = (^o 



1 - 


A 


2" 









; while asymmetrical, the ut — k relation in 

3 / , \ 41 



In addition, the embed- 



ded feature of the self-similarity in the RSJ model is also extended to obtain the power 
law relation between the threshold frequency of each branch, utu, and the branch number, 
n. The threshold frequency of each branch, lutu, is resca led by the branch number n as 



. According to the physical implication of the renormalization 
method conventionally, the exponent 5 has been found to present different kinds of behaviors 
under three stages of the damping constant A;, respectively. The generic relation connecting 
ujt and k of the above equations provide a crucial guide post to determine which kinds of 
solutions (simple or complicated) can be taken in accord with the strength of the frequency 
at a constant damping in the rf-biased nonlinear oscillators quantitatively. 
Finally, we do emphasize that the crucial formulas, ujt — ojq 
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Cl 



for potentials 

3 / , \4' 



possessing symmetrical property and ujt = cuq 
for asymmetrical potentials, can be tailored the following applications intuitively: (1) we can 
choose somewhat exact solutions of the system in a chaos-free region with a reference to the 
restriction of the parameters to obey such crucial formulas. For example, the symmetrical 
solution in the chaos-free region can be almost precisely expressed in the form of the single 
mode as X — a ■ sin(a;i -|- 0). (2) if we know how to exactly locate where the chaotic region 
is, then we can plan immediately concerning how to avoid the chaotic noise or to control 
the chaos in these nonlinear oscillating systems. 
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5. Figure Captions: 

Fig.l The state diagram of RSJ model with k = 1.0. The shaded area denotes the comph- 
cated behavior regions including those caused by various kinds of bifurcation. There 
is no chaotic behaviors found as the rf frequency exceeds ujt- The inset is the cor- 
responding potential well, 1 — cos x. The digits from 1 to 7 represent the number of 
the local wells which the motion runs. The phase portraits of the states marked from 
'a' to 'g' are presented in the figure 3. The symbols of "S" and "SB" represent the 
symmetrical and asymmetrical solutions, respectively. 

Fig. 2 There bifurcation diagrams of RSJ model with k = 1.0 and F from 0.0 to 8.0 and for 
cu — 0.313, 0.400, and 0.600, respectively. All the responses on the Poincare section are 
carried into the intervals of [— tt, tt] by the transform of x ± 2m7r.The similar scenarios 
of bifurcation appear itself again and again in each branch. 

Fig.3 The phase portraits with k = 1.0, oj = 0.313, and (a)F = 1.00, n = 1; (b)F = 1.75, n = 
2; (c)F = 2.60, n = 3; (d)F = 3.55, n = 4; (e)F = 4.50, n = 5; (f)F = 5.60, n = 6; 
(g)F = 6.50, n = 7; (h)u; = 0.05, F = 8.00, n = 51. The value of n represents the 
number of the local wells which the motion runs. 

Fig. 4a The numerical simulation data oi cut — k compared with the fitting function, cut — 
cuo[l — {k/koy], where cuq = 1.450 and ko = 1.681. 

Fig.4b The numerical simulation data cutu ~ n ioi various values of k and fitted by the 
functional relation cutu — <^c^~'^ ■ 

Fig. 4c The dependence on k of fitting constant cUc obtained from fig.4(b). 

Fig.4d The dependence on k of fitting constant 6 obtained from fig. 4(b). 
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Fig. 5 (a) The state diagram of asymmetrical one-well Duffing oscillator including the pri- 
mary (Aj) and subharmonic (Bj) resonance regions. C denotes the complicated 
solution region enclosed by period-doubling and reverse period-doubling transition 
boundaries and Es escaping boundary, cut indicated the maximum rf frequency 
for chaotic behaviors existing. The potential well is drawn in the inset, (b) The 
ujt — k relation, LUxik) = u!q[1 + ci{k/ko) — C2{k/koy + cs{k/ko)^ — C4{k/ko)'^] with 
cuo = 2.226, A;o = 2.200, ci = 0.101, C2 = 4.011, C3 = 4.495, and C4 = 1.585. The black 
circles denote the experimental simulation data and the solid line for fitting. 

Fig. 6 (a) The state diagram of symmetrical two-well Duffing oscillator including the primary 
{Aj) and subharmonic (Bj) resonance regions. C denotes the complicated solution 
region enclosed by period-doubling and reverse period-doubling transition boundaries 
and Hup the jump-up boundary. The potential well is presented in the inset. As long 
as the rf frequency is above Ut and the rf amplitude is below H^p transition boundary, 
there is no chaotic behaviors found, (b) The cut — k relation (9), lut — <^o[l ~ {k/ko)^] 
with luq — 3.218 and ko — 2.050, of the rf-biased Duffing oscillator with double well 
potential. 

Fig. 7 (a) The state diagram of symmetrical one- well Duffing oscillator including the primary 
(Aj) and subharmonic (Bj) resonance regions. C denotes the complicated solution 
region enclosed by period-doubling and reverse period-doubling transition boundaries 
and Eg escaping boundary. The potential well is presented in the inset. As long as the 
rf frequency is above lut, there is no chaotic behaviors found, (b) The lut — k relation 
(10), OUT — ouo[^ — [k/koY] with ouq — 1.528 and ko — 1.985, of the rf-biased Duffing 
oscillator with one well potential. 

Fig. 8 The state diagram of the Duffing oscillator with infinite potential well. The presented 
transition boundaries are including saddle- node (H), symmetry-breaking {AS), and 
period-doubling and reverse period-doubling (C) bifurcation. The shaded terrain de- 
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notes the complicated solution regions. There is no ujt found. 



16 



